subroutine init_fptr(nlog,nxec,nstp,be)
  use param
  use fptr
  implicit none
  integer,intent(in) :: nlog,nxec,nstp
  real(dbl),intent(in) :: be
  real(dbl) :: e,r,mine,maxe,td
  integer :: i,tn,pos
  ! Initialize xptr
  xptr=0
  read(nxec,rec=1) mine
  read(nxec,rec=2) maxe
  read(nxec,rec=3) td
  tn=int(td)
  lenx=2*tn+2
  do i=1,tn
     pos=2*i+2
     read(nxec,rec=pos) e
     if (be.le.e) then
        xptr=pos
        exit
     end if
  end do
  ! Initialize pptr
  pptr=0
  read(nstp,rec=1) mine
  read(nstp,rec=2) maxe
  read(nstp,rec=3) td
  tn=int(td)
  lenp=2*tn+2
  do i=1,tn
     pos=2*i+2
     read(nstp,rec=pos) e
     if (be.le.e) then
        pptr=pos
        exit
     end if
  end do
end subroutine init_fptr
